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Beams of polar molecules can be focused using an array of electrostatic lenses in alternating 
gradient (AG) configuration. They can also be accelerated or decelerated by applying an appropriate 
high voltage switching sequence to the lenses. AG focusing is applicable to molecules in both low- 
field and high-field-seeking states and is particularly well suited to the problem of decelerating heavy 
molecules and those in their ground rotational state. We describe the principles of AG deceleration 
and set out criteria to be followed in decelerator design, construction and operation. We calculate the 
longitudinal and transverse focusing properties of a decelerator, and exemplify this by 2D-imaging 
studies of a decelerated beam of metastable CO molecules. 
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I. INTRODUCTION 

During the last few years, a variety of techniques have 
been demonstrated to produce samples of trapped neu- 
tral molecules [l| . One of these exploits the force that a 
polar molecule experiences in an inhomogeneous electric 
field to change its motion. This force attracts molecules 
to regions of high or low electric field depending on the 
sign of the Stark shift. Some small polar molecules in 
low-field seeking states have been decelerated using a se- 
ries of pulsed electric fields. These include CO [2|, NH3 
and ND3 0, OH [1,[|, NH ih H2CO 0] and SO2 [^]- In 
the case of ND3 [J] and OH [5[ packets of Stark deceler- 
ated molecules have subsequently been electrostatically 
trapped. We aim to extend this deceleration method to 
heavy polar molecules including bio-molecules. Of par- 
ticular interest are molecules such as YbF which are be- 
ing used in experiments aimed at detecting time-reversal 
symmetry violating interactions leading to a permanent 
electric dipole moment (EDM) of the electron, which is 
a sensitive probe for physics beyond the Standard Model 
[3, [ifll- Decelerated molecules offer an increased sensi- 
tivity for these experiments. Stark deceleration of bio- 
molecules allows one to prepare samples of selected con- 
formers for further studies. 

Deceleration of heavy polar molecules is difficult for 
two reasons: (i) For a given velocity of the beam, the 
kinetic energy of the molecules, and thus the number of 
electric field stages required to bring the molecules to 
a standstill, is proportional to their mass, (ii) At the 
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electric field strengths required for deceleration, all low- 
lying rotational levels of heavy molecules have a negative 
Stark shift. In these states the molecules are attracted to 
a maximum of the electric field, i.e. to the electrodes. In 
order to guide high-field-seeking molecules through the 
decelerator dynamic focusing schemes need to be used 
which, typically, have an order of magnitude smaller ac- 
ceptance than the schemes used to guide low- field- seeking 
molecules. We have recently demonstrated a decelerator 
for high-field seeking molecules [ll|, [I^l using the alter- 
nating gradient principle. Those experiments showed the 
feasibility of the AG deceleration technique, which we 
discuss in detail here. 



Our paper is organized as follows. In Sec. |TI] we out- 
line the method of alternating gradient Stark decelera- 
tion, and set out some general principles that will guide 
us through the rest of the paper. In Sec. IIIII the Stark 
shift of polar molecules is discussed in more detail, draw- 
ing on the particular cases of metastable CO, YbF and 
benzonitrile. In Sec. IIVI we present three simple elec- 
trode geometries that may be used to make a single lens 
of the alternating gradient array and discuss the mer- 
its of these geometries. In Sec.|V]we consider the motion 
of the molecules through the decelerator, present the tra- 
jectories and phase-space distributions, and calculate the 
transmission of both idealised and real decelerators. In 
Sec. IVII we present an experimental study of the trans- 
verse focusing properties of an AG decelerator, by mea- 
suring the two-dimensional distributions of a decelerated 
beam of metastable CO molecules. Our results are com- 
pared to calculations. A summary of our main conclu- 
sions and a discussion of future prospects are given in 
SecEHl 
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ing one of these capacitors. Again, by switching the elec- 
tric fields at the appropriate times the molecule will be 
decelerated. A series of switched electric fields can thus 
be used to decelerate a pulsed molecular beam. The ini- 
tial velocity of a seeded supersonic beam of molecules is 
in the range 250-2000 m/s, depending on the mass and 
temperature of the carrier gas. Typically, more than 100 
stages are required to decelerate these beams to zero ve- 
locity. In order to have useful transmission, it is therefore 
of utmost importance that the trajectories through the 
decelerator are stable. 

For a force field, -F(r), to keep a particle in static equi- 
librium around r = 0, two conditions must be met. The 
applied force must vanish at r = 0, and, for small dis- 
placements, the force field should tend to restore the 
particle towards r = [531 ■ To achieve the latter it is 
necessary that the divergence of the force be negative, 
V • F < [5J|. The force acting on the molecules in an 



inhomogeneous electric field is given by: 



F{r) ^~VW{E), 



(1) 



with W{E) being the Stark shift of a polar molecule in 
an electric field of magnitude E = \E\. The properties 
of this force field were analyzed in a seminal paper by 
Auerbach, Bromberg and Wharton [1^. For molecules 
that experience a linear Stark shift in the applied field, 
W = —^effE, it was shown that: 



FIG. 1: (a) Layout of an alternating gradient decelerator for 
polar molecules showing the first four deceleration stages. 
Each electrode pair acts both to focus and decelerate the 
molecules, (b) Cross-section of a single lens formed from two 
20 mm long rods with hemispherical ends, 6 mm in diameter 
and spaced 2 mm apart. Potential energy along the z-axis for 
metastable CO molecules in the a"^n, J—l, f2=l, MQ—+1 
level, when the potential difference between the electrodes is 
20 kV. The high- voltage switching procedure is indicated: the 
voltages are turned on when the bunch of molecules reaches 
the 'HV on' position, and are turned off once they reach the 
'HV off position. 



II. 



GENERAL PRINCIPLES 



A Stark decelerator consists of a series of capacitor 
plates. A polar molecule that has its dipole oriented anti- 
parallel to the electric field (a low-field seeker) will gain 
potential energy when entering one of these capacitors 
and will therefore be decelerated. When leaving the field 
of the capacitor it will lose potential energy and so be ac- 
celerated back to its initial velocity. If we switch the field 
off before the molecule has left the capacitor, it will keep 
its lower velocity. Similarly, a molecule that has its dipole 
oriented parallel to the electric field (a high-field seeker) 
will accelerate when entering, and decelerate when leav- 



V-F = 



3 



t,j,k=l 



dxk 



9$\ 
dxij 



dxk 



^2$ 

dxidxj 
dxidx^ 



92$ 
dxkdxj 



(2) 



where $ is the electrostatic potential and ^eff is an ef- 
fective dipole moment which depends on the particular 
molecular state. Using Schwartz's inequality, it can be 
seen that the sum is always positive. Therefore, for 
molecules having a linear Stark shift the sign of V • F 
is determined solely by the sign of ^^ff [53 ■ Thus, for 
molecules that have a negative fXeff (low- field seekers), 
V • F < 0, and focusing is straightforward. For molecules 
that have a positive ^ieff (high-field seekers), V • F > 0, 
and focusing is more problematic. 

The difficulty of focusing high-field-seeking molecules 
is analogous to the situation for ions, for which V ■ F — 
qV • F = in free space, where q is the charge of the 
ion. Therefore, techniques routinely applied to ions can 
be translated to polar molecules. Three schemes are gen- 
erally employed, (i) Circular motion; In a cyclotron the 
curvature of the trajectory adds a force which, in an ap- 
propriately shaped magnetic or electric field, stabilizes 
the motion of the ions lj|. A similar stabilization can 



be achieved for polar molecules. For example, molecules 



that have a Unear Stark shift in an apphed field flying at 
a distance r from the axis of a capacitor formed by two 
coaxial cylinders, experience a force proportional to 1/r^. 
They therefore move in stable Kepler-type orbits around 
the central electrode. This technique has been used to 
focus molecules in high- field seeking states |15l - [l7| . (ii) 
Alternating gradient (AG) focusing] Alternating gradient 
focusing of charged particles was pioneered by Courant, 
Livingstone and Snyder [1^, ^^ and is now applied in vir- 
tually all particle accelerators. An AG array consists of a 
series of magnetic or electric quadrupole lenses that focus 
ions in one direction while defocusing them in the other 
direction. By alternating the orientation of these fields 
it is possible to obtain net focusing in both directions. 
As this stabilization is due to the motion of the ion it- 
self, it is referred to as 'dynamic' stability. Application of 
the technique to focus polar molecules was demonstrated 
experimentally by Kakati and Laine [20j|, by Giinther et 
al. [21-23] and by Bromberg [ll]. More recently, the AG 
technique was used to focus metastable argon atoms re- 
leased from a magneto-optical trap [2^ and cesium atoms 
in an atomic fountain |2q . Furthermore, the transmis- 
sion of methylfluoride molecules through a 15 m long AG 
beamline was modelled and optimized ^27J. (iii) Einzel 
lens; In an Einzel or uni-potential lens an ion is subjected 
to an acceleration along the axial direction followed by 
an equal deceleration. In the radial direction ions are fo- 
cused on entering the fringe field and defocused on leav- 
ing it. This results in a net (dynamic) focusing effect. 
The focusing is only effective when the change in kinetic 
energy of the ions is a substantial fraction of their ini- 
tial energy and so Einzel lenses are only useful for low 
energy ion beams. A similar effect is obtained for polar 
molecules entering and leaving a field region, but again 
its usefulness is restricted to very low energy beams, e.g. 
for loading slow molecules into a trap. 

For a decelerator (or accelerator) for polar molecules 
alternating gradient focusing seems the obvious choice. 
Fig. [TJa) shows the general form of the experimental 
setup. The AG lenses are formed from a pair of cylin- 
drical electrodes to which a voltage difference is applied. 
Molecules will be defocused in the plane containing the 
electrodes while being focused in the orthogonal plane. 
As the molecules move down the beamline the focus- 
ing and defocusing directions alternate. The defocusing 
lenses have a smaller effect on the molecules than the 
focusing lenses, not because their power is smaller (it 
is not), but because the molecules tend to be close to 
the axis when they are inside the defocusing lenses and 
further away from the axis when they encounter the fo- 
cusing lenses. Molecules in high-field seeking states are 
accelerated while entering the field of an AG lens and are 
decelerated while leaving the field. By simply switching 
the lenses on and off at the appropriate times, AG focus- 
ing and deceleration of polar molecules can be achieved 
simultaneously. Fig. [Ijb) shows the potential energy (the 
Stark shift) along the z-axis of a single lens for a represen- 
tative high-field seeking molecule. The molecules enter 



each lens with the electric fields turned off so that their 
speed is unchanged as they enter. The fields are then 
suddenly turned on, and the high-field seeking molecules 
are decelerated as they leave the lens and move from a 
region of high field to one of low field. This process is 
repeated until the molecules reach the desired speed. As 
indicated in the figure, the amount of deceleration can 
be controlled by choosing how far up the potential hill 
the bunch of molecules has climbed before the fields are 
turned off (i.e. by moving the 'HV off' point in Fig.[TJb)). 
Similarly, the effective length of each lens, L^g, can be 
controlled by varying the amount of time that the fields 
are on (i.e. by moving the 'HV on' point in Fig. [TJb)). 

A prototype machine of this type has been used to 
decelerate high-field seeking metastable CO molecules 
from 275 m/s to 260 m/s JTj. More recently, an im- 
proved device decelerated ground state YbF molecules 
from 287 m/s to 277 m/s, corresponding to a 7% reduc- 
tion of the kinetic energy [ij]. Since then, YbF and 
CaF molecules have been decelerated using an array of 
21 lenses at Imperial College London, while at the Fritz- 
Haber-Institut in Berlin, CO and benzonitrile molecules 
have been decelerated using an array of 27 lenses. These 
results will be presented elsewhere. 



III. THE STARK SHIFT IN POLAR 
MOLECULES 

We now discuss the Stark shift in more detail for 
some representative polar molecules. The Stark shift of 
a molecule is a function of the electric field magnitude 
E = \E\. It is useful to define a dimensionless parame- 
ter, A, that describes the strength of the electric field. In 
the context of an idealized rigid-rotor molecule [23| with 
dipole moment /x and rotational constant B (in energy 
units), the appropriate dimensionless ratio is A = fiE/ B. 
When the electric field is 'weak', A ^ 1, the Stark shift is 
quadratic in A and can be calculated using second order 
perturbation theory. The states are best labelled by the 
rotational angular momentum quantum number J, and 
its projection, M, onto the field axis. If A is increased, 
states of different J are increasingly strongly mixed un- 
til, in the stroiig-field limit (A ^ 1), the states are called 
'pendular' |29l-l31]. In that case they are labelled by the 
quantum numbers Vp and M, with Vp = 2 J — \M\, and 
states of the same Vp but different M are degenerate. In 
this strong- field limit, all the low- lying states are high- 
field seekers. This limit is of most interest for our present 
discussion. Within this high-field, pendular state model, 
the Stark shift, W, is given by 



Wivp, \)/B = -A + {vp + l)(2A)i/2, 



(3) 



showing that the Stark shift becomes asymptotically lin- 
ear in the electric field. 

We find it useful to define the effective dipole moment, 
/iejj, more generally as 
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FIG. 2: Stark shifts and effective dipole moments for (a) tlie 
J=l, fl= 1 levels of the a^II state of CO, (b) the lowest four 
rotational levels in the ground state of YbF and (c) the lowest 
eight rotational levels of benzonitrile. In (a) and (b) the basis 
set includes all levels up to J— 10, while in (c) all levels up to 
J— 30 were included. 
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which converges with increasing electric field to its maxi- 
mum value /i, reached when the body-fixed dipole mo- 
ment is parallel to the external electric field. From 
Eq. ^ and (jj]) the effective dipole moment is given, in 
the strong-field limit, by 



IJ-eff =M 



(5) 



Note that the effective dipole moments of states hav- 
ing the same values of Vp converge once the strong-field 
criterion (A ^ 1) is met. By contrast, the convergence of 
all the ^eff to a single value is a much slower one, scaling 
as X~^/^ and so requiring A^/^ » 1 for low values of Vp. 
In strong fields the effective dipole moment varies little 
with applied field, and for small changes of the field it 
can be approximated as a constant. This is a very useful 
approximation in the context of an alternating gradient 
lens where the field is high and does not vary greatly 
across the aperture of the lens. 

Figure [SJa) shows the Stark splitting and the effective 
dipole moments for the J=l, ^2= 1 level of the a'^II ex- 
cited state of CO. The electronic ground state of CO has 
a small dipole moment (O.lDebye) and rotational levels 
in this state only experience a second order Stark effect in 
realizable fields. By contrast, the metastable a^II state of 
CO (lifetime 3.7 ms) has a dipole moment of 1.37Debye 
(IDebye is equivalent to 0.0168 cm^^/kV/cm). Being a 
H-state, all the rotational levels are doubly degenerate. 
As the separation of the nuclear motion and the elec- 
tronic motion is not exact, this degeneracy is lifted and 
each rotational level is split in zero electric field into two 
levels with opposite parity. For the J=l, ^= 1 level this 
A-doublet splitting is A=394MHz. The two A-doublet 
levels are coupled by an electric field, leading to levels 
with a mixed parity that have non-zero space-fixed elec- 
tric dipole moment. The Stark shift of the two A-doublet 
levels in a small static electric field of magnitude E is 
found by diagonalizing the energy within a single rota- 
tional manifold: 
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where J denotes the total angular momentum, while fl 
and M are the projections of J onto the body fixed and 
space fixed axes, respectively. At higher electric fields, 
the Stark effect includes coupling to states of the same M 
but different J. Since low-lying states have many states 
of higher J above them, this coupling ultimately turns 
them all into high-field seekers. For example, the upper- 
most level in Fig. [H^a) is weak-field seeking, but at fields 



above ^-^400 kV/cm it becomes high- field seeking. The ef- 
fective dipole moments for the lowest rotational levels of 
CO are shown below the Stark-curves in Fig. ^a). Cal- 
culation of the Stark shift in metastable CO is discussed 
in detail by Jongma et al. [S^l ■ 

Figure [^b) shows the energy and effective dipole mo- 
ments of YbF in the X'^T,'^ electronic ground state as 
a function of the electric field strength [33]. The states 
are labelled by the rotational quantum number N and its 
projection onto the electric field axis, Mn- This Stark 
effect is caused by the mixing of rotational levels. The 
ground rotational state is high-field seeking at all fields. 
Other states, such as the N = 1, Mm = state are low- 
field seeking at small electric fields but become high-field 
seeking at larger field values. For the N — l,Mjv — 
state, the turning point occurs at an electric field of 
~ BB/n, corresponding to only ~18kV/cm for the heavy 
YbF molecule. For the fields in Fig. [D^b), the high-field 
condition A 3> 1 is satisfied, and one sees that the ^eff 
values converge for states of the same Vp = 2J— \M\, but 
different M. It is also evident in the figure that conver- 
gence of the effective dipole moments to the single value, 
fi, is very slow, as discussed earlier. 

In Fig. [2][c) we show the Stark effect and effective 
dipole moments in the lowest rotational states of benzoni- 
trile, calculated using experimentally determined con- 
stants '34l|. Benzonitrile is an asymmetric top, and there- 
fore levels with the same J are mixed by the electric field 
as well as those having AJ = ±1. As a molecule of this 
size has rather small rotational constants, all rotational 
levels become high-field seeking in relatively weak elec- 
tric fields. The jumps between dipole moment curves in 
Fig. Hfc) are caused by avoided crossings. Details on the 
calculation of the Stark shifted energy levels in an asym- 
metric top molecule can be found elsewhere J35| . 

Table U gives the relevant properties for Stark decel- 
eration for a selection of polar molecules. These prop- 
erties are the Stark shift, effective dipole moment, rota- 
tional constants and mass. Values are given at a field 
of lOOkV/cm and are for molecules in the electronic and 
rovibronic ground state (with the exception of metastable 
CO). The number of electric field stages required to bring 
molecules with a certain initial velocity to rest depends 
on the ratio of their Stark shift to their mass. The focus- 
ing properties of molecules fiying at a certain velocity de- 
pend on the ratio of their effective dipole moment to their 
mass. One can see from the table that these molecules, 
though widely different in mass, have similar ratios of 
effective dipole moment to mass, and of Stark shift to 
mass, and so will be focussed and decelerated similarly. 
The Stark shifts and effective dipole moments are gen- 
erally dependent on the specific quantum state that the 
molecule is in. For many experiments one would like to 
decelerate molecules in a variety of quantum states simul- 
taneously. This can be done when all the rotational states 
have the same dependence on electric field, as is more or 
less the case for the polyatomic aromatic molecules with 
small rotational constants listed in the table. For com- 



plex molecules such as tryptophan, the decelerator offers 
the intriguing possibility of selecting a specific confor- 
mational isomer out of the various conformers known to 
co-exist in a supersonic beam J4l|, as the individual con- 
formers have distinctly different values oi fieff- 



IV. ELECTRODE GEOMETRY 
A. The field of an infinitely long lens 

Although the electrodes of a decelerator are in short 
segments along the beam direction, the basic focusing 
properties are best elucidated by first considering the 
case of long electrodes. In this section we discuss how 
to design a set of electrodes that minimizes the aberra- 
tions of an AG lens. In an aberration-free lens, molecules 
experience a harmonic interaction potential in the trans- 
verse plane. As discussed in Sec. [TTl the potential will fo- 
cus along one direction and defocus along the other. For 
molecules that experience a linear Stark shift the ideal 
form for the field strength is also harmonic; E{x,y) — 
Eq + r]{x'^ ~ y^)- As we shall see, this field cannot be re- 
alized but it is possible to produce a field that is a good 
approximation to this ideal one. We follow a similar ap- 
proach to that given in j27| . 

In a region devoid of charges the electric field can be 
derived from the electrostatic potential (^ as E — — V<I>, 
with V^$ = 0. In 2D, $ may be represented by a multi- 
pole expansion as: 



$(x,2/) = $0 
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Here r = ^/{a 



y 



'■) and 9 — tan "'^ (f ) are the usual 



cylindrical coordinates. a„ and 5„ are dimensionless con- 
stants. r-Q and $o are scaling factors that characterize the 
size of the electrode structure and the applied voltages, 
respectively. The electric field ma gnitude at the centre 
is given by Eq = ($o/'''o)-\/ai + ^i- The n = 1 terms in 
Eq. ([7]) represent a constant electric field, while the n ^ 2 
and n = 3 terms represent the familiar quadrupole and 
hexapole fields that have been used extensively to focus 
molecules in low-field seeking states 
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Equation ^ represents the most general form of the 
electrostatic potential consistent with Laplace's equation. 
Now we choose the coefficients to be suitable for mak- 
ing a good lens. We require the magnitude of the elec- 
tric field to be non-zero at the origin, and symmetric 
under reflection in the x- and y-axes. To achieve this, 
we make <& symmetric under reflection in the a;-axis and 
anti-symmetric under reflection in the y-axis by setting 
all bn = and retaining only the terms of odd n. Antic- 
ipating the result that high-order terms only introduce 
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Rotational state 


Stark shift 


Effective dipole 
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Mass 








(cm-l) 
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at 100 kV/cm 


at 100 kV/cm 


A/ B / C 




CO (a^ni) [32j 




\j = i,Mn = -i) 


-1.25 


0.0135 


- / 1.68 / - 


28 


CaF [36, 37] 




\J = l/2,Am = +l/4> 


-3.43 


0.0420 


- / 0.34 / - 


59 


YbF [33] 




\J = l/2,AfQ = +l/4> 


-4.91 


0.0569 


- / 0.24 / - 


193 


ND3 [38] 




\J = 1,MK = -1) 


-1.27 


0.0134 


- / 5.14 / 3.12 


20 


pyridazine [39| 




\Jk^kJM\) = \OooO) 


-5.59 


0.0624 


0.21 / 0.20 / 0.10 


80 


benzonitrile [34] 




\Jk^kJM\) = \OooO) 


-6.71 


0.0711 


0.19 / 0.051 / 0.040 


103 


tryptophan f40| 


I 


\Jk^kJM\) = \OooO) 


-6.25 


0.0646 


0.041 / 0.013 / 0.012 


216 




II 




-4.72 


0.0494 


0.039 / 0.014 / 0.012 






III 




-1.71 


0.0183 


0.033 / 0.017 / 0.013 






IV 




-11.68 


0.120 


0.032 / 0.016 / 0.013 






V 




-12.28 


0.126 


0.043 / 0.011 / 0.0096 






VI 




-11.37 


0.116 


0.045 / 0.011 / 0.0095 





TABLE I: A selection of polar molecules with their relevant properties for AG focusing and deceleration. 



undesirable non-linearities into the force we choose to re- 
tain only ai, a^ and a^. Hence: 



$(a;, J/) = $0 I ai— + as^ ^ 
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(8) 



From this potential, we obtain the electric field magni- 



tude, E{x,y) = ^(|f )' -f (||) . Throughout the re- 

gion r < To this can be expanded as a power series in 03 
and 05. For the case 05 <C 03 ^ oi we obtain 
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The first two terms have the desired form and domi- 
nate the expansion. The other terms produce focus- 
ing aberrations. It might appear advantageous to set 
os/ai = (03/01)2/3 so as to cancel the cross term. How- 
ever, the best policy is to minimize a^/ai [27i | because 
the x^ + y^ term is also a damaging aberration. 

To produce these fields we need to choose electrodes 
whose surfaces map onto the equipotentials. We are free 
to choose 03/01 either positive or negative and will dis- 
cuss both cases in turn. For example. Fig. ^a) shows 
equipotentials for 03/01 = +1/7 and 05 = 0. The choice 
of 03 is constrained by the condition 03 <C Oi while re- 
maining large enough to provide significant focusing. The 






x/ro 



FIG. 3: Electrostatic equipotentials of equation ((S]) for the 
cases (a) as/ai = -1-1/7, as = and (b) aa/ai = —1/7, 
as = 0. Red and blue colourings correspond to positive and 
negative potentials respectively. The white lines show elec- 
trode surfaces designed to follow these contours. 



solid white lines show electrodes, chosen to be circular 
for ease of construction, that match the equipotentials 
closely. They have radii of i? = Srp and are centred at 
X = ±4ro, leaving a gap of 2ro. Because these elec- 
trodes do not match the equipotential exactly, higher 
order terms appear in the field. From a fit to the nu- 
merically calculated electrostatic potential we find for 
this geometry 03/ai = 0.143 and 05/03 = 0.143. It is 
noted that this two-rod field can be solved analytically 
and that 03/01 = (ro/i?) / (2 + ro/R) and a^/a^ — 03/01, 
in agreement with our fit. 

The two charged rods are schematically depicted in 
Fig. mja), while Fig. Hljb) shows the magnitude of the 
electric field they produce as a function of distance along 
the X-axis (solid line) and y-axis (dashed line). A high- 
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FIG. 4: The case of a two-rod lens, (a) Schematic view of 
the lens formed by two 6 mm diameter rods spaced 2 mm 
apart with a potential difference of 20 kV. (b) The electric 
field strength versus displacement along the a;-ajds (solid line) 
and the t/-axis (dashed line), (c) Forces on a CO molecule in 
the high-field-seeking component of the a^II, J=l, n=l level. 
Dashed line: focusing force, Fy{0,y). Solid line: defocusing 
force Fx{x, 0). The sign of this force has been reversed for ease 
of comparison. The gradient of both lines near the origin is 



fc=0.37cm ^/mm^ 



field seeker will be defocussed along x and focussed along 
y. Fig. \^c) shows the corresponding forces exerted on 
CO molecules in the high field seeking component of the 
a'^n, J—1, n,—l level. The dashed line is the restoring 
force generated by a displacement along y, and is seen to 
be roughly linear. The solid line indicates the defocusing 
force along x. The sign of this force has been reversed 
so that the two can be compared directly. The two have 
equal gradients near the origin. Further away from the 
origin, the non-linearity due to 05 acts to strengthen the 
defocusing power whereas the focusing is weakened. We 
will see later that this difference reduces the acceptance 
of an alternating gradient decelerator. 

The rather large value of 05/03 in this two-rod case can 
be reduced by adding two grounded electrodes tangential 
to the $ = equipotentials at y = ±2.65ro. These are 
shown dashed in Fig. [3]^ a), where for simplicity we have 
given the new rods the same radius R. In this case the 
coefficients become a^/ai = 0.157 and 05/03 = 0.070. At 
the expense of a slightly less ideal field one can position 
the four identical electrodes at the corners of a square. 
This has the advantage that one is free to choose in which 
plane the field focuses or defocuses by simply switching 
the voltages [4^ . Using electrodes of radius ro with their 
centers placed on the corner of a square of side 3ro yields 
o-z/o,! = 0.59 and a^/a^ = 0.056. The rather large value 
of 03/01 introduces higher order terms in the field, even 
though 05/03 is quite small. A disadvantage of this field 
geometry is that the electric field strength on the beam 
axis is only half that on the electrodes. This makes the 
configuration less suited for use in a decelerator as the 
energy removed per stage is proportional to the central 
field. The geometry is useful for guiding molecules, as 
was recently demonstrated by Junglen et al. [44| . 

We turn now to the case of negative 03/01 illustrated 
in Fig. [Sjb) where we have chosen 03/01 = —1/7 and 
05 = 0. This is well approximated by electrodes of radius 
R = 2.3ro, with a minimum gap of 2ro as shown by the 
solid white lines of Fig. Eljb). The precise field produced 
by these electrodes has 03/01 — —0.139 and a^/a^, — 
—0.014. This geometry compares very favourably to the 
cases considered in Fig. [DJa) with regard to minimizing 
05 and hence the lens aberrations. As this geometry is 
symmetric, we can again reverse the focus and defocus 
directions very easily by interchanging the potentials on 
the top-right and bottom-left electrodes of Fig.[21[b). The 
field at the centre is 41% of the maximum, which is dis- 
advantageous for a decelerator. This geometry was dis- 
cussed by Liibbert et al. [23] in the context of focusing 
ICl molecules. 



B. End effects 

Until now we have assumed the electrodes to be in- 
finitely long. We now discuss the infiuence of end effects. 
Eq. ([2]) gives the divergence of the force on a molecule 
with linear Stark shift. If we restrict our attention to the 
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FIG. 5: The 'force constants' for CO (a^D, J=l, ^=1, 
MQ,=+1) along the molecular beam axis and near the exit 
of the lens shown in Fig. [Hb) and Fig. [IJa). Solid line: de- 
focussing constant k^ = —dF^/dx. Dashed line: focussing 
constant ky = —dFy/dy. Dash-dotted line: axial gradient 
k^ = -dF^/dz. The three are linked by Eq. ^^. The gray 
shaded area indicates the region of the lens's hemispherical 
end. 



axis of our beamline, this equation is greatly simplified. 
Consider, for example, the pair of electrodes shown in 
Fig. m The electrostatic potential is symmetric about 
the a:-z-plane and the y-z-plane is one of antisymmetry 
with $ = everywhere on this plane. All the electrode 
geometries considered here have this property. On the 
beamline, the intersection of these two planes, it follows 
that all the derivatives in equation ([2]) are zero apart 
from d^/dx and d^^jdzdx. Since —d^/dx is the only 
non-zero electric field component, its magnitude is the 
total electric field strength E. Thus, on axis, 



V-F = 



E 



dE 
dz 



(10) 



Inside the lens dE/dz = and dFz/dz = and it follows 
that the spring constants in the two transverse directions 
kx = —dFx/dx, ky = —dFy/dy are equal and opposite, 
as shown in Fig. |4ljc). The equality fc^; = —ky means 
that the focusing and defocusing powers are equal. In 
the fringe field of the lens dE/dz ^ and dF^/dz ^ 
and we find that 



kx + ky 




-^^eff 



^d ( IdE 
dz \E dz 



d^E 
dz"^ 



(11) 



Due to the inhomogeneity of the electric field along z , the 
defocusing force becomes larger than the focusing force 
near the exit of the lens, whereas the focusing force is 
larger than the defocusing force further away from the 
lens. This is illustrated in Fig. [SJ that shows the three 



force constants kx,ky and kz for the high- field seeking 
CO molecules as they approach the end of a lens formed 
by two rods with hemispherical ends. In this figure, the 
origin of z is at the point of inflection {d'^E/dz'^=0) and 
the region of the hemispherical ends is indicated by the 
grey shaded area. We begin with the left hand side of the 
figure where the end-effects are negligible. Here, kz = 
and therefore k^ = —ky = 0.37cm~^/mm^ the same as 
in Fig. m As the molecules approach the exit of the lens, 
they experience a decelerating force which can be seen 
in the figure as a positive kz- According to Eq. ([TU]), 
this is accompanied by a corresponding decrease in the 
sum kx + ky. We see this as a strengthening of the defo- 
cusing constant kx near the exit, which becomes nearly 
twice as strong as the focusing constant. As we will see 
in Sec. IV Al this is another mechanism, in addition to 
the aberrations discussed in Sec. IIVAI that significantly 
reduces the transmission of an AG decelerator. A more 
gradual termination of the rods, e.g. a prolate spheroid 
replacing the hemisphere, reduces both the first and sec- 
ond derivatives of E with respect to z and so reduces 
end-effects in accordance with Eq. ((TT|) . In comparison 
with the two-rod configuration, the four-rod arrangement 
of Fig. ini^b) is also found to have a more favourable field 
at the exit of the lens. 



V. MOTION OF THE MOLECULES THROUGH 
THE DECELERATOR 

In this section we investigate the motion of molecules 
through the decelerator and discuss criteria for optimiz- 
ing the transmission. The first part of this section deals 
with transverse stability, and the second part with longi- 
tudinal stability. This division is based on the assump- 
tion that the transverse and longitudinal motions can be 
treated independently. This is an approximation whose 
validity we discuss at the end of the section. 



A. Transverse motion 

In describing the transverse motion, we start by assum- 
ing that the molecules experience a linear force that fo- 
cuses them along one direction and defocuses them along 
the other. The orientation of successive lenses alternates. 
The lenses have lengths L and are separated by drift re- 
gions of length S where the molecules experience no force. 
For molecules moving with a constant velocity, Vz , along 
the molecular beam axis, the equation of motion in a lens 
can be written as d'^x/dz^^K^x = 0, where the plus sign 
applies in a focusing lens, and the minus sign applies in 
a defocusing lens. The number of oscillations per unit 
length inside a focusing lens is K/27r and is related to 
the force constant khy n — \J\k\/'mv1. We also define 
the angular oscillation frequency i7 which, for the linear 
Stark effect reads 




(12) 



For a molecule with initial position x{zo) and velocity 
v{zo), the equation of motion can be written as 



from our use of {x,Vx) as state variables, rather than 

{x,Vx/Vz). 

The Courant-Snyder parameters and the phase ad- 
vance are related to one another: 



a(z) 



Vz dfi{z) 
' 2 dz ' 



(17a) 



x{z) 

Vx{z) 



= M{z\zo) 



x{zo) 
Vx{za) 



(13) 



The transfer matrix Af (zJzq) is then given by 



M(z|zo) = 



cos K,l 
— Q. sin k1 
1 l/v, 
1 
cosh K.I 
51 sinh nl 



cos kI 



i sinh kI 
cosh kI 



F: 



O: 



D: 



7(z) 



l + a^{z) 



$ = — 



z+l. 



focusing lens 
drift space 
defocusing lens 



P{z'\ 



■dz'. 



(17b) 



(17c) 



where I = z — zq- The transfer matrix is written as F in 
a focusing lens, as £) in a defocusing lens and as O in a 
drift region. 

The transfer matrix for any interval made up of subin- 
tervals is just the product of the transfer matrices of the 
subintervals: 



Equation (J17b|) ensures that the matrix has unity deter- 
^Minant. When expressed in this form, the transfer ma- 
trix acquires an extremely useful property, namely that 
the matrix describing N lattice units is identical to the 
matrix for a single unit, but with $ replaced by N^. 
Equation (I17cp shows that $ is independent of z, since 
the integral is taken over one complete period of the peri- 
odic function /3. Note that the stability criterion becomes 
— 1 < cos$ < +1 and so is satisfied if $ is real. 

The trajectory of a molecule moving through the ideal 
lattice is given by 



M{z2\zo) = M{z2\zi)M{zi\zo). 



(15) 



x{z) = y//3{z)eiCos{(l){z) + 6i), 



(18) 



A single repeating unit of the alternating gradient ar- 
ray has the transfer matrix F{L).OiS).D{L).0{S). We 
have written the lengths L and S explicitly here, but will 
usually drop them. The transfer matrix for an array of 
N such units is M = {FODO)^ . Ahernatively, it can 
be useful to introduce frequent deceleration sections into 
longer lenses using a configuration M = {FO)"{DO)". 
This structure with n = 3 is used in the decelerator that 
we present in Sec. |Vll In order for molecules to have 
stable trajectories through any such array it is neces- 
sary that all the elements of the transfer matrix remain 
bounded when N increases indefinitely. This is the case 
when -1 < ^Tr{M) < +1 (see, for example, fist). 

It is useful to parameterize the transfer matrix of one 
repetitive unit with length Iceii as [l9[ 



M{z + lcen\z)^ 



cos ^ + a sin $ 



(sin$ 



-7sin$ cos$ — asin$ 



(16) 



where a{z), /3{z) and 7(2:) are z-dependent parameters 
with periodicity equal to that of the lattice and are 
known as the Courant-Snyder parameters. $ is known 
as the phase-advance per cell. Note that I3{z) and 7(2:) 
are expressed in seconds and 1/seconds, respectively, 
rather then in meters and 1 /meters as is customary in 
the charged particle accelerator literature. This follows 



where e^ and Si define the initial conditions of this partic- 
ular molecule, and ^(z) is a z-dependent phase given by 
(j){z) — 1/vz L l//3(z') dz'. Equation ([T5|) shows that the 
motion is a product of two periodic functions, the first of 
wavelength Iceii and the second of wavelength 2'Klceii/'^- 
When <I> ^ 27r, the first motion has a short wavelength 
and is known as the micromotion, while the second has a 
much longer wavelength and is called the macromotion. 
This motion is identical to that of an ion in an rf trap 
[4a |. From x{z) and the relationships that hold between 
the Courant-Snyder parameters, it can be shown that 



7(z)a;^ -I- 2a{z)xVx + (3{z)v'^ — 



(19) 



This equation defines an ellipse in the phase-space 
whose coordinates are x and Vx ■ The shape of the ellipse 
evolves periodically with z, but always has the same area 
TT£i. A set of molecules having many different values of 
6i but the same value of e^ will all lie on the same ellipse. 
Furthermore, a distribution of molecules with all possible 
values of Si and all values of e^ in the range < e^ < e, 
will all lie inside the ellipse characterized by e. Again, the 
shape of this ellipse evolves periodically, but its area is a 
constant, ire. The value of e defines the size of the beam 
in phase-space, and is called the emittance of the beam. 
Equation P^ tells us that the transverse displacements 
of a set of molecules lie within a beam envelope given by 
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FIG. 6: (a) Trajectories of metastable CO molecules flying 
with a forward speed of 315 m/s through a (FO)^ {DO)^ array. 
The parameters are K=38.7m~^, L=2mm and 5'=28mm, 
corresponding to $ = n/6. The shaded area shows the 
beam envelope bounded by ±-^/^, for a constant aperture 
of d=2mm. (b) As (a), but with L=6mm and 5'=24mm cor- 
responding to $ = 7r/2. (c) Phase-space area occupied by 
the beam at four positions in the unit cell with the position 
spread in mm and the velocity in m/s. (i) Centre of the defo- 
cusing set, (ii) between the defocusing and focusing sets, (iii) 
centre of the focusing set and (iv) between the focusing and 
defocusing sets. 



the periodic function ±A//3(z)e. The velocity spread hes 
within a beam envelope given by zLy^j{z)e. 

In Fig. mja), some trajectories are plotted for 
metastable CO molecules travelling at 315 m/s through 
an alternating gradient array of type {FO)^{DO)^. The 
parameters of the array are k = 38.7m~^, L = 2mm 
and S — 28 mm. Here, the phase-advance is 7r/6, and 
the micromotion of wavelength Iceii is superimposed on 
a macromotion whose wavelength is 12 1 ceii- The shaded 
area of the figure shows the envelope of the transmit- 
ted molecular beam as it passes through the array. For 
such small values of the phase-advance, there is only a 
small difference between the maximum and minimum 
sizes of the beam envelope. As the phase-advance in- 
creases, the modulation of the beam envelope increases. 
This is demonstrated in Fig. Hfb) which shows trajecto- 
ries and beam envelope for the same value of k but with L 
increased to 6 mm and S decreased to 24 mm. The phase 
advance is now 7r/2, meaning that molecules return to 
their starting point after Alceii- Figure [Sljc) shows the 
phase-space distribution of the beam at four positions 
within the unit cell. In graph (i), the molecules are at 
the centre of the defocusing triplet. Here, the transverse 
size of the beam is at its minimum. The beam is diverging 
as it enters the focusing lenses (ii) , and reaches its maxi- 
mum size at the centre of the focusing triplet (iii). Graph 
(iv) shows that the beam is converging when it enters the 
defocusing lenses. The fact that the transverse size of the 
beam is larger in the focusing lenses than in the defocus- 
ing lenses, and that the forces are proportional to the 
off-axis displacements, accounts for the stability of the 
array. Since the stability relies upon the motion itself, it 
is commonly referred to as 'dynamic' stability. 

We now calculate the transverse acceptance of the AG 
array. If we suppose that the array of lenses has a uniform 
aperture d throughout its length, then the beam whose 
emittancc is e will be transmitted without loss provided 
that the envelope fits inside the aperture, i.e. provided 
that yj I3{z)e < d/2 everywhere in the array. The trans- 
verse acceptance is the phase-space area occupied by the 
beam of largest emittance consistent with this criterion. 
This area is nd"^ /{4:(3max)- From Fig. [6] we see that /? is 
always a maximum at the centre of a focusing lens. To 
calculate the transverse acceptance of a lattice, we sim- 
ply find the value of /3 at this position using Eq. p^ 
with M = F{L/2) ■ 0{S) ■ D{L) ■ 0{S) ■ F{L/2). From 
this we find that /3 is a dimensionless number divided by 
il and so the acceptance is a multiple oi d^Vt. d and fi 
are the natural scaling parameters of the problem, and 
trajectories are invariant when plotted in (x / d, Vx / (dfl)) 
space. 

In Fig.[7I^a) and Fig.[71[b) we show the transverse accep- 
tance calculated for lattices of FODO and {FOf{DOf 
cells, respectively. The acceptance is plotted as a func- 
tion of the two dimensionless parameters that define the 
lattice, kL and kS. The acceptance (in either transverse 
direction) is given in units of d^Q,. One sees that the 
highest transverse acceptance is O.lSGd^ri and is obtained 
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FIG. 7: Acceptance in one transverse direction of an infinitely 
long alternating gradient array, as a function of kL and nS. 
The acceptance is in units of d Q. Also shown are some con- 
tours of $ which define the stability region, (a) A FODO 
array, (b) A (FO)^{DO)''^ array. Straight solid and dashed 
lines, cross, triangle and dot are all referred to in the text. 



when kL ~ 1 and S <Si L. By contrast, a single, infinitely 
long focusing lens has an acceptance of (7r/4)d^f2, which 
is over four times larger. The figure also shows some con- 
tours of $, filling the region of stability bounded by the 
cos<i> = ±1 contours. The cos$ — +1 contour corre- 
sponds to the vertical line k,L = 0. 

In the experiments discussed in Sec. IVIl the physical 
structure is fixed aX L + S = 30 mm, but the effective 
length of the lenses can be varied by adjusting the high- 
voltage switch-on time (see Fig.[TJb)). Thus, the possible 
operating conditions lie on a straight line. Two such lines 
are shown in Fig. [71[b), indicating the operating condi- 
tions for the experiments where metastable CO was used 
at forward speeds of 630 m/s and 315 m/s. The cross 
and triangle placed on Fig. El^b) correspond to the set- 
tings used to calculate the trajectories in Fig. El^a) and 



Fig. EJb), with small and large phase advance respec- 
tively. 

It is worth considering the scaling behaviour of the 
transverse acceptance with aperture, d. The accep- 
tance along each direction scales as d^Jl, and J7 scales 
as EJ d~^ (Eq. ([H])). It is natural to operate the decel- 
erator at the maximum field that can be achieved, which 
is determined by the breakdown field. If the transverse 
scale of the lenses is increased, with Eq held constant by 
corresponding increases in the applied voltages, the ac- 
ceptance along each direction is linear in d, implying that 
one should make the aperture as large as possible. There 
is, however, a practical upper bound Vmax on the applied 
voltages. Once this value of Vmax is reached, Eq scales 
as d^^ and the acceptance in each transverse direction 
is proportional to d^/"^ . Although the acceptance contin- 
ues to increase with d, the decreasing value of Eq results 
in an undesirable decrease in the energy loss per decel- 
eration stage. Furthermore, increases in the acceptance 
cease to be useful once the transverse emittance of the 
beam is fully contained within the transverse acceptance 
of the decelerator. 

We have seen that transverse stability depends on 
maintaining coherence between the oscillation of the 
molecules and the structure of the array. This makes 
the alternating gradient focussing particularly sensitive 
to deviations from the ideal, such as nonlinear terms in 
the force, end effects, and misalignments of the lens ar- 
ray. We now discuss the impact of each of these on the 
transverse acceptance. Fig. [5] shows the two-dimensional 
transverse acceptance of an {FOY{DOY array calcu- 
lated for the various lens geometries discussed in Sec.lIVI 
The lens-lens spacing, L -I- 5, is fixed at 30 mm, and the 
operating conditions correspond, for metastable CO, to 
the two straight, solid lines drawn on Fig. El^b). Figure 
[8l[a) is calculated for CO molecules with a forward speed 
of 315 m/s, while Fig. |Sl[b) is for a speed of 630 m/s. 
The calculation uses the known Stark shift of metastable 
CO i32,] and the electric fields obtained from SIMION 
[431 . The trajectories of, typically, 510^ molecules with 
random initial positions and velocities are traced through 
a 96-lens array by numerical integration. Line (i) shows 
the acceptance obtained for a set of perfect linear lenses 
with Eq = 95kV/cm, a^/ai — 1/7 and tq = 1mm. The 
lens aperture, d, is taken to be 2.1mm for reasons that 
will become clear shortly. Since this line is for perfect 
lenses, it sets the scale for the forthcoming cases. Line (ii) 
shows the acceptance obtained when the force has a small 
non-linearity, corresponding to the real field of the four 
lens geometry of Fig. [3]Ja) , with the two high voltage elec- 
trodes held at ±10 kV and the other two grounded. Once 
again, ro is 1mm and this gives i?o=95kV/cm. End- 
effects are not considered. The curve displays structure 
that is absent in the ideal case, and the region of high 
acceptance is seen to be narrower. We chose the value of 
d in the ideal case (curve (i)) so as to give the same maxi- 
mum acceptance. This suggests the definition of an effec- 
tive aperture, deff, 2.1 mm in this case. It is interesting to 
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FIG. 8: The acceptance along the sohd hnes indicated in 
Fig. [7] for (a) CO seeded in Xenon (uz=315ni/s) and (b) 
CO seeded in Argon (wz =630 m/s). Lenses with four differ- 
ent electrode geometries have been used, as indicated in the 
legend. For comparison, the transmission of perfectly linear 
lenses with a minimum aperture of 2.1 mm is also shown. 



note that this effective aperture is shghtly larger than the 
real 2 mm gap between the two electrodes at high volt- 
age. This occurs because the beam envelope is smaller in 
the defocusing direction than in the focusing direction, 
and the smaller gap is in the defocusing direction. The 
2D transverse acceptance scales as d^^ . We next increase 
the size of the non-linear contributions to the force, by 
removing the two grounded electrodes of Fig. |31Ja). As 
discussed in Sec. IIVI this approximately doubles the ra- 
tio 05/03. All other parameters are kept constant, and 
end-effects are not yet considered. The acceptance in this 
case is given by line (iii) . We find the impact of the non- 
linearities to be very detrimental indeed. The effective 
aperture is reduced to 1.6mm. As shown by line (iv), 
a further reduction in acceptance occurs when we intro- 
duce the fringe-field aberrations at the entrance and exit 
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FIG. 9: The calculated acceptance for an array of 96 ideal 
lenses, as a function of the degree of misalignment (defined 
in the text). Random (i) and constant (ii) misalignments of 
a FODO array with K=50m~^, L=S=15mm. Random (iii) 
and constant (iv) misalignments of a {FO)'^{DO)'^ array with 
K=38.7m^^, L=6mm and S=24mm. 



of each lens. Here, the electrodes of the two-rod lens 
have hemispherical ends of radius 3 mm as outlined in 
Sec. lIVBl The effective aperture for this case is 1.3 mm. 
We have also considered the four-electrode geometry of 
Fig-E^b), with the electrodes at ±10 kV and r^ = 1mm, 
giving Eq = 47kV/cm. As discussed in Sec. IIVI this 
geometry results in very small non-linearities, the ratio 
05/03 being a factor of 10 smaller than in the two-rod 
case. The calculated acceptance is shown by line (v) and 
is divided by a factor 4 for ease of comparison. As Eq, 
and hence $7 is smaller in this case, the maximum accep- 
tance is shifted to higher values of L. The effective aper- 
ture for this configuration is 2.9 mm making it by far the 
most effective configuration considered. Note that this 
effectiveness is due to the small value of 05/03 and not 
due to the sign of 05/03; indeed, calculations show that if 
we reverse the sign the acceptance is the same. As men- 
tioned in Sec. IIVI a disadvantage of this geometry is that 
the field on axis is only half that of the two-rod geome- 
try, for the same maximum electric field. A decelerator 
composed of these 4-rod lenses would therefore require 
twice as many deceleration stages as one composed of 2- 
rod lenses. The much improved acceptance comes at the 
cost of increased decelerator length. 

Finally in this section, we discuss how the inevitable 
misalignments of a real machine affect the transverse ac- 
ceptance. We consider two types of misalignment. In 
the first type, referred to as 'random', the centre of each 
lens is displaced horizontally and vertically from the axis 
by amounts chosen at random from normal distributions 
with full- width at half-maximum Ar. In the second type, 
referred to as 'constant', lenses that focus in the hori- 
zontal direction are perfectly aligned to one another and 
define the axis, while the lenses that focus in the vertical 
direction are all displaced from the axis by an amount 
Ar in both transverse directions. These types of mis- 
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alignment tend to occur naturally in the construction of 
the decelerator. For example, in our decelerators, each 
electrode is mounted into one of four common bars to 
which the high voltages are applied. The degree of ran- 
dom misalignment is determined by the machining pre- 
cision and the construction technique. Misalignment of 
the four bars relative to each other results in a constant 
displacement of the horizontal lenses from the vertical 
lenses. 

We have calculated how the transverse acceptance of 
the decelerator diminishes as the degree of misalignment 
increases. Sensitivity to misalignments is found to de- 
pend somewhat on the values of k, L and S. Some rep- 
resentative cases are shown in Fig. [HI This figure gives 
the 2D transverse acceptance for an alternating gradi- 
ent array of 96 ideal lenses (linear force, no end effects) 
as a function of the degree of misalignment, Ar/r^, r^ 
being the radius of the circular aperture defined by the 
lenses. Each line in the figure has been normalized to 
the acceptance obtained for perfect alignment. Line (i) 
gives the result in the case of random misalignments in a 
FODO array with K=50m^^, L=15mm and 5=15 mm. 
One sees that random misalignments of Ar ~ O.OSra are 
sufficient to reduce the acceptance by 50%. Line (ii) plots 
the effect of a constant misalignment for the same param- 
eters, showing this to be a less severe misalignment in this 
case. A 50% drop in acceptance is reached when this mis- 
alignment reaches Ar ~ 0.05ra. Lines (iii) and (iv) plot 
the random and constant cases for a {FOY{DOY array 
with K=38.7m^^, L—Qmra and S'=24mm. Here, the ac- 
ceptance is not so sensitive to the random misalignments, 
while the sensitivity to the constant misalignment is the 
same as for the FODO case. These curves give some typ- 
ical scenarios. In general, we find that the transmission 
is less sensitive to misalignments for smaller values of kL, 
and that misalignments are most severe when the array 
is operated close to the stability boundary at $ = tt. 

In our decelerators with ra= 1mm, we have achieved 
values of approximately 20 /im for the size of the ran- 
dom misalignment by specifying tight machining toler- 
ances where appropriate. An alignment jig was used to 
reduce the constant type of misalignment below ~50/im. 



electric fields when the molecules have not yet left the 
field of a lens completely, as shown in Fig. [TJb). This 
ensures that molecules at the head of the pulse lose more 
kinetic energy, while those at the tail lose less. In this 
way, molecules with a suitably narrow spread of longitu- 
dinal position and velocity can be be confined to a small 
area of phase-space throughout the decelerator. This be- 
haviour, known as phase stability, has been discussed ex- 
tensively in the context of decelerating weak-field seeking 
molecules [48]. 

To analyze the longitudinal motion, we begin by Taylor 
expanding the on-axis potential energy of a single decel- 
eration stage (see Fig. [TU)) around the point of inflection 
at z = 0: 



W{z) = W{Q) + W'{Q)z + W"'{Q)z^ 



"OF" 



(20) 



and we have 



where W'{Q) ^ ^l^^g, W"'{Q) 

used the fact that d'^W/dz'^\z=o = 0. Close to z = the 
potential energy can be approximated using the first two 
terms only. The switching sequence is constructed such 
that a hypothetical molecule, the so-called 'synchronous 
molecule', always reaches the same position Zs of the rel- 
evant lens at the moment when the fields are turned off. 
The change in kinetic energy of the synchronous molecule 
is the same in every lens, Wions — W{zs), W\cns being the 
Stark shift of the molecule inside the lens. 

When the energy taken out per stage is small com- 
pared to the total kinetic energy of the molecules t48| . 
i.e., when Au ^ v, one can describe this change in en- 
ergy as originating from a constant force 



F, = 



Wi, 



W(Zs) 



(L + S) 



(21) 



where {L + S) is the distance that the synchronous 
molecules travels between two subsequent switching 
times. The difference of the force on a non-synchronous 
molecule at position z and the force on the synchronous 
molecule at position Zg is now given by 



B. Longitudinal motion 

In order to decelerate or accelerate the molecules, time- 
varying electric fields are applied. A molecule in a high- 
field seeking state will gain kinetic energy as it enters the 
field of a lens, while it loses kinetic energy as it leaves the 
lens, as shown schematically in Fig. [ijb). If the electric 
field is switched on while the molecule is inside a lens 
there is no change to its kinetic energy but the molecule 
will decelerate as it leaves the lens. The moment when 
the field is switched on determines the effective length 
Leff of the lens and hence the focusing properties. The 
moment when the field is switched off determines the 
deceleration properties of the lens. We switch off the 



F-F, 



-W{z) + W{zs) 
(L + S) 



W 



(L + S) 



(z-zs). (22) 



Consequently, the non-synchronous molecules oscillates 
around the synchronous molecule with an angular fre- 
quency given by 



W 



m{L + S) ' 



(23) 



with m being the mass of the molecules. 

In our experiment with metastable CO molecules, 
described in Sec. I VII L + S'=30mm, and 
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FIG. 10: The Stark shift of metastable CO (a^Di, w'=0, 
J'=l, Mn=l) molecules as a function of their position along 
the molecular beam axis. The dotted lines indicate the posi- 
tion of the synchronous molecule when the electric fields are 
switched on and off for two different settings of the decelerator 
(expressed as the energy change per stage, AVK, compared to 
the maximum possible energy change, ISWmax)- In the lower 
part of the figure some (closed) trajectories in phase-space are 
shown for non-synchronous molecules, relative to the position 
and velocity of the synchronous molecule. 



M^'(0)=0.29cni~^/mm resulting in a longitudinal 
frequency a;2/27r=320Hz. 

Figure [TU] shows the trajectories of a few non- 
synchronous molecules, plotted in phase space relative to 
the position and velocity of the synchronous molecule. 
From these numerical simulations, non synchronous 
molecules are found to oscillate around the synchronous 
molecule with a frequency of a;z/27r=330Hz, close to the 
frequency given by Eq. (P5|) . The thick curves in Fig. [TU] 
show the outermost trajectories of molecules that are still 
phase stably decelerated. The longitudinal acceptance is 
about 50mm-m/s, when ISW — Q.llSWmaxi and three 
times larger when IS.W — ^flhlSWmax- 

Alternating gradient deceleration can also be applied 
to low-field seeking molecules. In this case, a slightly 
more complicated switching pattern must be used to 
achieve both longitudinal and transverse stability. Sup- 
pose we want the synchronous molecule to lose an en- 
ergy ^.hlSWmax per stage. The fields should be turned 
on well before the synchronous molecule approaches the 
lens, and turned off again when it is half way up the po- 
tential hill. This ensures that non-synchronous molecules 



oscillate around the synchronous one as before. Once the 
molecules are well inside the lens, the fields are turned 
on again to focus the molecules, and must be turned off 
before they approach the exit to ensure that they are not 
accelerated out of the lens. The first high voltage pulse 
determines the amount of deceleration whilst the second 
determines the effective length of the lens. 



C. Coupling between the longitudinal and 
transverse motion 

The transverse stability depends on the longitudinal 
velocity because k is inversely proportional to v^,. Sup- 
pose the experimental settings at the start of the de- 
celerator correspond to the point (kL, kS) indicated by 
the dot in the transverse acceptance plot of Fig. [71 As 
the molecules are decelerated, their position on this plot 
moves away from the origin along the dashed line. Even- 
tually, this point will move out of the region of stabil- 
ity, and the beam will be lost. To avoid this, either 51 
or L and S must be altered along the array in sympa- 
thy with the decreasing speed. A decrease in J7 could 
be achieved by decreasing the curvature of the electric 
field. This could be done without altering the on-axis 
field, which governs the energy loss per stage. However, 
the transverse acceptance is proportional to fi, and un- 
less the decrease in S7 is compensated by an increase in d, 
this will lead to beam loss. A more satisfactory approach 
is to decrease L and S so that Ljv^ and Sjv^ remain 
constant. With this approach the lenses will be long at 
the beginning of the decelerator and since deceleration 
occurs only at the end of each lens the overall length 
may then become undesirably large. That problem can 
be circumvented by splitting each lens into several parts, 
i.e., by replacing the FODO array with the more general 
{FO)"^{DO)" array. As the velocity is decreased, n is 
also decreased until, at the end of the decelerator, n — 1. 
In this way the beam can be decelerated stably to a small 
fraction of its initial speed. 

In the longitudinal direction, molecules oscillate 
around the position of the synchronous molecule, caus- 
ing the effective lens length experienced by a molecule 
to vary according to the phase of its longitudinal oscil- 
lation. This also couples the longitudinal and transverse 
motions, possibly leading to parametric amplification of 
the transverse oscillation |49|]. The coupling can be sup- 
pressed by designing the lenses to be long compared to 
the longitudinal spread of the decelerated beam. 



VI. 2D IMAGING OF AN AG DECELERATED 
BEAM OF CO MOLECULES 

In order to demonstrate the performance of the alter- 
nating gradient decelerator, we have carried out exper- 
iments on CO molecules in the a'^ll state. The main 
reason for choosing metastable CO molecules for these 
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FIG. 11: Scheme of the experimental setup. CO molecules are prepared by UV laser excitation to a high-field seeking level 
of the metastable a'^II state and pass through an array of 12 lenses arranged in the {FO)^{DO)^ configuration, with the last 
group of electrodes oriented vertically. The transverse distribution of high-field seeking metastable CO molecules is recorded 
20 cm after the decelerator using a micro channel plate, phosphor screen and CCD camera. 



experiments is that (i) they can be prepared in a single 
quantum state at a well-defined position and time, and 
(ii) their velocity distribution as well as their transverse 
distribution can be readily recorded. A scheme of the 
experimental setup is shown in Fig. 1111 A pulsed beam 
of CO is produced by expanding a mixture of CO with 
either Xe or Ar into vacuum, using a modified solenoid 
valve. When seeded in Argon the mean velocity of the 
beam is i'i=630m/s, corresponding to an initial CO ki- 
netic energy of Efcj„ = 480 cm^^. When Xenon is used 
and when the valve housing is cooled to 180 K (just above 
the boiling point of Xe at the pressure used) the mean 
velocity of the CO molecules in the beam is reduced to 
315 m/s (mixture of 20% CO in Xe) or to 275 m/s (mix- 
ture of 5% CO in Xe) corresponding to initial kinetic en- 
ergies of 120 cm^^ and 89 cm^^ , respectively. In all cases, 
the velocity spread is approximately 10%, corresponding 
to a translational temperature of about 1 K. 

The metastable CO molecules are prepared in a sin- 
gle quantum state by direct laser excitation on the spin- 
forbidden a^H (u'=0) <— X^T,^ (^v"—0) transition, using 
narrow-band pulsed 206 nm (6.0 eV) radiation. In the ex- 
periments reported here, the laser is tuned to excite the 
lower A-doublet component of the J'—l a^Hi level via 
the i?2(0) transition. By setting the polarization of the 
laser perpendicular to the stray electric fields present in 
the excitation region only the MH. = 1 high-field seeking 
level is prepared. 

The CO molecules pass through a 1.0 mm diame- 
ter skimmer into a second, differentially pumped, vac- 
uum chamber housing the 35 cm long AG decelerator. 
The decelerator consists of 12 equidistant 20 mm long 
lenses, separated by 10 mm long drift regions. The 
lenses are arranged in four groups of three, with the 
first group of electrodes oriented horizontally and the 



last group of electrodes oriented vertically. The lenses 
are formed from two circular electrodes (bold white lines 
of Fig. [31Ja)), with ro=lmm and i?=3mm, and have 
hemispherical ends. The two opposing rods are simul- 
taneously switched between OkV and ±10 kV by two in- 
dependent high voltage switches. The electric field on 
the axis is 95kV/cm, corresponding to a Stark shift of 
-1.2 cm~^ for the metastable CO molecules. The Stark 
shift on the molecular beam axis is shown as a function 
of z in Fig. [TJb). 

The molecules land on a micro-channel plate (MCP) 
detector placed on the beam axis. The 6 eV energy of 
the excited state is sufficient to release Auger electrons 
from the surface. These are amplified and detected on a 
phosphor screen using a CCD camera (LaVision GmbH). 
Thus, the 2D distribution of the metastable CO beam is 
recorded. The detection efficiency of the MCP detector 
is estimated to be about 10~^ [50|. Detection efficiencies 
> 10% can be obtained by letting the molecules impinge 
on a flat gold surface kept at 500 K and redirecting the 
Auger electrons towards an MCP mounted off-axis. Un- 
fortunately, our attempts to build sufficiently distortion- 
free optics to image the electrons from the gold surface 
onto the MCP detector failed. Therefore, the longitudi- 
nal characteristics of the decelerator are recorded using 
the gold plate detector, but the 2D distribution of the 
decelerated CO is measured with the MCP directly inter- 
cepting the beam. The initial intensity of the metastable 
CO beam is monitored simultaneously by detecting the 
a'^n <— X^E^ fluorescence near the entrance of the AG 
decelerator with a photomultiplier tube (PMT). 

Figure [12] shows the measured time-of-flight (TOF) 
distributions for several values of the energy change per 
stage, AW. The timing sequences were chosen so that 
the decelerator would always act on a group of molecules 
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FIG. 12: Observed time-of-flight (TOF) distributions of 
metastable CO molecules over the 54 cm path length through 
the apparatus, for 4 different settings of the decelerator (ex- 
pressed as the energy change per stage, AVF). The lowest 
curve shows the TOF-distribution when the electric fields are 
off. The measurements (thick curves) have been given an 
offset for clarity. The thin curves show the results of 3D tra- 
jectory calculations. The vertical dashed line indicates the 
expected arrival time of a molecule flying with a constant 
velocity of 275 m/s. 



with initial speeds centred on Vi — 275 m/s. The lowest 
curve is the TOF distribution obtained with no voltages 
applied to the decelerator. Molecules with a speed of 
275 m/s arrive at the time indicated by the dashed line in 
the figure. Using the 12 stages, the speed can be reduced 
to 263 m/s, or increased to 286 m/s, depending on the 
timing of the switched fields. In these experiments the de- 
celerator was operated at ±8 kV. The thin curves showing 
the results of the trajectory calculations discussed earlier 
describe the TOF distributions for the decelerated bunch 
rather well. Similarly good agreement is obtained in the 
simulation of the accelerated beam (not shown). The ra- 
tio of the time-integrated signal with decelerator on and 
off is ^1.5 as also predicted by our simualtions. This ra- 
tio is about 10 times larger than measured in a previous 
experiment on metastable CO |11| . 

Figure [13] shows transverse distributions of the 
molecules, measured 20 cm downstream from the decel- 
erator exit for a variety of effective lens lengths. The 
imaging detector is switched on for a short period (20 fis 
when Ui=630m/s and 40 /iS when Vi=315m/s) so that 
only the decelerated molecules are detected. The images 



are formed from ^10'* detected molecules accumulated 
over many shots (2000 when Vi— 630 m/s and 10^ when 
Vi—315 m/s). The false colour is a measure of the number 
of molecules detected in each pixel. A calculated spatial 
distribution is shown beside each experimental image. 

We discuss first the data with Vi = 630 m/s, Fig. lT^ a). 
The top image shows the 2D distribution obtained 
when there are no voltages on the decelerator. As 
expected, this measurement simply shows equally dis- 
tributed molecules within the aperture formed by the 
electrodes. The detected distribution has a width-to- 
height ratio of about 1.3 because the last 3 electrode pairs 
are oriented vertically and end 9 cm closer to the detec- 
tor than the last horizontal electrode pairs. By looking 
at the outline of the electrodes, one can also see that the 
lens array was not perfectly aligned. Indeed, the centre 
of the horizontal lens set is displaced by 250 ^m from the 
center of the vertical set (we have since reduced this mis- 
alignment to less than 50 fim) . This misalignment has 
been included in the simulations. 

The remaining images in Fig. I13f a) show how the pro- 
file of the beam changes as the effective length of the 
lenses is increased. We consider the vertical direction 
first. In this direction, the last set of lenses are defo- 
cussing, and the beam is diverging when it exits the de- 
celerator. One might expect this divergence to increase 
with increasing lens length. However, the smaller size 
of the beam inside the dcfocussing lenses (compare (a) 
and (b) in Fig. [5]) tends to compensate for the increased 
power of those lenses. As a result, the vertical divergence 
changes very little in our experiments and the height of 
the distribution is approximately constant in all the im- 
ages. In the horizontal direction, the last set of lenses 
are focusing lenses, and so the beam is converging when 
it exits the decelerator. In contrast to the defocussing 
lenses, the size of the beam inside the focussing lenses 
is fixed, being determined by the lens aperture. As the 
lens length increases the exiting beam converges more 
strongly and so the width of the detected distribution 
decreases. When the lens length is about 8 mm, a focus 
is formed in the plane of the detector. Here, the settings 
correspond to the dot placed on the transverse accep- 
tance plot of Fig. [7{b), where the phase-advance is 7r/6. 
The focus that is formed is rather aberrant, resembling 
a cross rather than a vertical line. As we will see shortly, 
the aberration is caused by non-linearities in the trans- 
verse forces. When the length is increased beyond 8 mm 
the focus lies in frontupstream of the detector and the 
horizontal width begins to increase again. In all cases, 
the experimental images agree very well indeed with the 
simulations. 

Turning now to the data with Wi=315m/s (Fig. ll3r b)). 
we see that the lower forward speed results in the focus 
being formed in the plane of the detector for a smaller 
value of Leff. Using a thin lens approximation, the fo- 
cal length is I/k^L. Since halving Vi, doubles k we ex- 
pect the focus to be formed for Leif=2mm, exactly as 
observed in the data. These settings correspond to the 



17 



(a) 



CO in Ar(V|=630m/s) 

Expenment Simulation 



(b) COinXe(V|=315m/s) 

Experiment Simulation 



10 



12 



14 





10 mm 



> <■ 



10 mm 



10 mm 



> <■ 



10 mm 



FIG. 13: Observed and calculated 2D distributions of a decelerated beam of CO molecules with an initial speed of (a) 630 m/s 
and (b) 315 m/s, for various values of the effective lens length (in mm). 
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(c) including end effects 




FIG. 14: Simulated 2D distributions 20 cm downstream from 
the decelerator exit, for CO molecules with Vi = 630 m/s. (a) 
Linear force, (b) true force, excluding end effects, and (c) true 
force. 



cross placed on Fig. [Tjb) , where again, the phase advance 
is 7r/6. As Leff is increased beyond 2 mm, the horizon- 
tal focus lies progressively further in front of the detec- 
tor and so the detected distribution grows progressively 
wider. When Lgjj > 10 mm, the trajectories become un- 
stable and there is a very sudden drop in the intensity 
of the simulated data. The experimental images show a 
less rapid drop in intensity, CO molecules being observed 
near the molecular beam axis, even when L^g = 14 mm. 
In addition, a diagonal cross shape is observed in the 
experimental data. This we attribute to molecules so 
strongly focused in the first three lenses that they es- 
cape and subsequently fly outside the decelerator. Our 
simulation could not follow such trajectories because the 
field used was bounded 4 mm from the beam axis and 
molecules outside this area were considered lost. We sur- 
mise that these molecules also give rise to the intensity 
observed at the centre of the images. 

By integrating the intensity of the measured distribu- 
tions one obtains the acceptance as a function of the ef- 
fective length of the lenses, analogous to the calculations 
shown in Fig. [S] However, a frequency drift of the UV 
excitation laser during the measurements compounded 
by the rather low counting rates, resulted in integrated 
beam fiuctuations of more than 50%, preventing a clear 
comparison of measured and calculated acceptances. 

Figure [U] reveals the origin of the aberrations present 
at the focus. We repeated the calculation for the case of 
Vi = 630 m/s, Leff = 8mm. In (a), the calculation was 
performed assuming a perfectly linear force, and we find 
the focus to be a vertical line without aberration. The 
aberrations appear when we include the non-linearities 
in the force (Fig. [TiT b)). showing this to be the primary 
factor in degrading the image quality. The addition of the 
end-effects degrades the image quality a little further, as 
shown in (c). 



VII. SUMMARY AND CONCLUSIONS 

In this paper we have presented the principles of al- 
ternating gradient deceleration of polar molecules along 
with the criteria that govern decelerator design and op- 



eration. We began by showing that heavy molecules in 
low-lying rotational states seek strong field, and so can- 
not be confined by static electric fields. Alternating gra- 
dient focusing can be used to stabilize the trajectories of 
such molecules in a Stark decelerator. We showed how 
the electric fields required to achieve alternating gradi- 
ent focusing can be produced using simple electrode ge- 
ometries. The simplest of these, the two-rod geometry 
shown in Fig. IDJa), has already been used to decelerate 
metastable CO and YbF molecules [lllli^, and was used 
in the imaging experiments presented here. This elec- 
trode geometry gives rise to significant non-linear terms 
in the force. These can be reduced, without reducing 
the on-axis field, by the addition of a pair of grounded 
electrodes. The non-linearities are further reduced when 
four high voltage electrodes are arranged as in Fig. [H^b), 
but at the cost of reducing the field on the axis. The 
fringe-fields of the lenses, essential for deceleration, tend 
to increase the defocusing power relative to the focus- 
ing power near the entrance and exit of the lens. This 
detrimental effect can be reduced by terminating the rods 
more gradually. The four-rod geometries are also more 
favourable than the two- rod geometry in this respect. 

We discussed the motion of the molecules through the 
decelerator in detail. When the transverse forces are 
perfectly linear, the transverse motion is best described 
using the formalism first set out in the context of the 
alternating gradient synchrotron [19|, and outlined in 
Sec. IV Al The trajectories are described in terms of the 
phase-advance and the envelope function. The transverse 
phase-space distribution is an ellipse whose shape evolves 
periodically through the lattice, but whose area is a con- 
stant. Figureiniillustrates the evolution of the beam enve- 
lope and the phase-space ellipse. We calculated the trans- 
verse acceptance of an array of ideal lenses as a function 
of kL and nS (Fig. [7]) and found the maximum accep- 
tance to be 0.186(i^ri, where d is the transverse aperture 
and 17 given in terms of the mass, the effective dipole mo- 
ment and the electric field curvature by Eq. ([T^ . Non- 
linearities in the transverse forces, and the fringe-fields 
of the lenses, severely reduce the transverse acceptance 
(Fig. [5]). Furthermore the degree of lens-to-lens align- 
ment required to achieve a high transmission was found 
to be demanding but attainable with high-precision ma- 
chining. 

The longitudinal motion was discussed in terms of a 
simple model of phase-stability. The most significant 
coupling between transverse and longitudinal motion is 
that the transverse motion depends on the forward speed. 
As the speed is reduced the molecules spend more time 
in each lens. A good way to handle this is to use an 
(FO)"(DO)" structure with n larger at the entrance 
than at the exit of the decelerator. 

We studied the focusing properties of an alternating 
gradient decelerator experimentally using an array of 12 
lenses, by measuring 2D images of a decelerated beam of 
metastable CO molecules. Trajectory simulations of this 
experiment reproduce the experimental findings. Non- 
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linearities in the force and effects due to rounded ends of 
the electrodes need to be included to obtain close agree- 
ment. These experimentally verified simulations predict 
a transverse acceptance of 2 (mm.m/s)^, see Fig. [51 and 
a longitudinal acceptance of 50mm.m/s, see Fig. [TOl 
These numbers can be compared to those that have been 
presented earlier for light molecules in low-field seeking 
states. In the Stark deceleration of ammonia, for in- 
stance, a transverse acceptance of 160 (mm.m/s)^ and a 
longitudinal acceptance of 10 mm.m/s has been obtained 
Q. As discussed above, the transverse acceptance is very 
much smaller than for ideal lenses and we can expect 
a more sophisticated lens design to yield a tenfold im- 
provement. A further increase of the acceptance could be 
achieved by increasing the transverse aperture, although 
a corresponding increase in the applied voltages would 
be needed in order to maintain the same on-axis electric 
field. 

In discussing the alternating gradient decelerator, we 
have also laid out the principles of a guide for high-field- 
seeking molecules. Unlike a decelerator, a guide does 
not need to be divided into segments along the beamline, 
and so can be free of end-effects. The symmetric four 
electrode geometry of Fig. ^h) makes an ideal guide be- 
cause its aberrations are small, its acceptance is high, and 
the focus and dcfocus directions are very easily switched. 
Another interesting application is to use this geometry 
as an m/ fj,eff filter, the equivalent of an m/q filter for 
ions [46] . The resolution of such a filter can be increased, 
at the expense of the acceptance, by tuning the focus- 
defocus duty-cycle away from 50%. In one plane, the de- 



focusing lenses are then longer than the focusing lenses 
making the stability region narrower and so increasing 
the resolution. Calculations indicate that a resolution 
A{m/ fieff)/i'm/ fieff) ^ 0.1 can be obtained at the cost 
of a factor of 4 in acceptance relative to the maximum. 

The decrease in velocity achieved in an AG decelerator 
so far has been rather small, but since the trajectories of 
the molecules through the decelerator are inherently sta- 
ble, no additional losses are expected when the number 
of stages is further increased. For the molecules listed in 
Table U and for many others, approximately 100 electric 
field stages are sufficient to bring them to rest. These 
molecules could subsequently be stored in a storage ring 
[sH or in an AC-trap |53 |. 
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